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Star cluster survival and compressive tides in 
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ABSTRACT 

Gravitational tides are widely understood to strip and destroy galactic substructures. 
In the course of a galaxy merger, however, transient totally compressive tides may 
develop and prevent star forming regions from dissolving, after they condensed to 
form clusters of stars. We study the statistics of such compressive modes in an N- 
body model of the galaxy merger NGC 4038/39 (the Antennae) and show that ~ 15% 
of the disc material undergoes compressive tides at pericentre. The spatial distribution 
of observed young clusters in the overlap and nuclear regions of the Antennae matches 
surprisingly well the location of compressive tides obtained from simulation data. 
Furthermore, the statistics of time intervals spent by individual particles embedded in 
a compressive tide yields a log-normal distribution of characteristic time r ~ 10 Myr, 
comparable to star cluster formation timescales. We argue that this generic process 
is operative in galaxy mergers at all redshifts and possibly enhances the formation 
of star clusters. We show with a model calculation that this process will prevent the 
dissolution of a star cluster during the formation phase, even for a star formation 
efficiency e as low as ~ 10%. The transient nature of compressive tides implies that 
clusters may dissolve rapidly once the tidal field switches to the usual disruptive mode. 
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1 INTRODUCTION 

Galaxy interactions are known to trigger star cluster for- 
mation. A typical example is NGC 4038/39 (the An- 
tennae) , where some 100 star clusters h as been de- 
tected (|Mengel et all l2005l : IWhitmore et al.1 120071 '). How- 
ever, it is unclear how many of them will survive in- 
ternal destruction; for instance, low star formation ef 
ficiency (SFE) or strong stellar feedback could d e- 
stroy the molecular cloud (|Elmegreen fc Efremovl [l997l) 
and lead to "in f ant-mortality" llBastian fc GoodwinI 200' 
IWhitmore et al.1 l2007l : Ide Griis fc Parmentieij l2007h . Fur 
thermore, these questions are central to interpreting the 
enhanced specific frequenc y and metall icity bi-modality of 
clusters in giant ellipticals (|Harridll99il ). In this Letter, we 
stress that totally compressive tides widely exist in galaxy 
mergers and slowdown the early dissolution of clusters. 

HST V-band and deep IR images of the Antennae galax- 
ies have suggested a near power-law cluste r mass function 
(^[M] (X M-^ for IC^ Mq < M < 10^ M^i lFall et a"i]|2005l : 
lAnders et all |2007| ) and an age distribution with a rough 
power-law fit (dN/dr oc t~^) suggesting that many young 



clusters dissolve rapidly after they are born inside their pri- 
mordial host gas cloud (Fall et al. 2005). 

The evolution of (broad-band) mass-to-light ratios of- 
fers a strong hint that young clusters are severely per- 
turbed and may ev e n be out of dynamical equi l ibrium 
(|Bastian fc GoodwinI l2006l : Ide Griis fc Parmentieij |2007| ). 
Yet the low inter-cluster velocity dispersion measured in 
the Antennae galaxy merger indicates that cluster formation 
proceeds from a coherent coalescence of large-scale material 
as opposed to in situ s hock dissipation and fragmentation 
jWhitmore et al.l l2005l ). Computer simulations of galaxy 
mergers have long suggested that streams of material (made 
up of both gas and stars) developing during a merger will 
form dense struc tures, so increasing the star- and c l uster- 
formation rates dBarne^ Il988l: iBarnes fc HernauistI Il992l : 



iMihos fc Hernauistlll99ilNaab et al.ll2006l ). A Schmidt law 
or variation thereof can be applied to circumvent numerical 
resolution issues and derive estimates of the star formation 
rate in rough agreement with obser vations (see e.g. Barned 
I2OO4I for an application to the Mice; Idi Matteo et al 1120071 1? 

Despite this wealth of progress, the long-term fate of 
star forming regions (bound structures or not) remains 
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largely unclear, especially with regard to star cluster com- 
plexes such as those of the Antennae galaxies. The situation 
is exacerbated in the case of a galaxy merger because the 
gravitational tides vary rapidly in time along with the mor- 
phology of the system as a whole. An observational snapshot 
showing a candidate star-forming aggregate lacks informa- 
tion about the duration and motion in space to determine 
whether the aggregate will remain bound or not. When mat- 
ter (gas or stars) is distributed in space in fragmented and ir- 
regular fashion, as is typical of galaxy mergers, the tidal field 
may exert either a disruptive outward stretch on an extended 
mass; or, on the contrary, it may exert an inward compres- 
sive force, which then cocoons the volume enclosed and in- 
creases its binding energy. While converging streams have 
long been sus pected of driving the forma t ion of tidal dwarf 
galaxies (e.g. 'Barnes fc HernquistI Il992l : iBournaud fc Dri3 
[2006; We tzstc in ct al. 2007), the impact of compressive tides 
themselves on sub-galactic scales has never been consid- 
ered strictly as a gravitational effect. Previous authors anal- 
ysed cluster formati on as the result of a gl obal trigger, us- 
ing sticky-particles (|Bournaud et al.| '2008') or SPIf simu- 
lations with a th reshold in pres sure dBckki fc Chiba 2002; 
iBekki et al]|2002l ') or in density l|Li et all 12003 ) for the for- 
mation of stars. 

Here, we take a different stand by exploring the im- 
pact of compressive tides on cluster-size volume elements. 
We present in §2 a high-resolution numerical model of the 
Antennae galaxy merger to match kinematic and optical 
data. In §3, we discuss the key results obtained from the 
computer simulation. We close with a discussion of how the 
process outlined here bears on the statistics of observed clus- 
ter populations. 



2 MERGER SIMULATION 

We picked two equal- mass bulge/disc/halo galaxy models to 
simulate an A ntennae-like major merger similar to that of 
iBarned ([1988") . We set up self-gravitating models using the 
magalie soft ware included with the NemcQ stellar dynami- 
cal package (|Teubenlll995l . version 3.2.4) . The equation s of 
motion were integrated with gyrfalcON (|Dehnenll2002l ) us- 
ing a Plummer softening kernel with h = 0.01 in units where 
the gravitational constant G = 1, a unit of time = 25 Myr, a 
mass unit = 3.6 x 10^" Mq and a length unit of 4.4 kpc. The 
total mass of each galaxy is equal to 7 mass units, where 
the exponential disc scale- length and mass rd = Md = 1. 
The bulge is a lHern guist (1990) model with a scale param- 
eter a — 1 and total masi|3 Mt — 1. Thus our progenitor 
galaxies are "SO - Sa" galaxies with a Toomre parameter 
Q — 1.5 which ensures local stability to axially symmet- 
ric fragmentation modes. An isothermal dark matter halo 
of total mass Mh = 5, of equal truncatio n radius and scal e 
length Th.t = Th = 7, was set up as in iHernauistI l|l993h . 
The mass ratio between individual component particles of 
^ 2 keeps two-body heating low. A total of 700, 000 particles 
were used for each galaxy, for a mass resolution of 3.50, 1.75 

^ http:/ /bima. astro. umd.edu/nemo/ 

^ Our m odel galaxie s differ mainly from Barnes's in that he used 
a lighter [Kind lll966[) model for the bulge. 



and 4.38 x 10^ Mq respectively for each componen10. Lower 
resolution simulations with 175,000 and 350,000 bodies per 
galaxies did not change significantly the statistics reported 
here, demonstrating that the results are robust to Poisson 
noise. No gas was included in our simulation. The typically 
small gas mass fraction would not affect the gravitation- 
nal potential significantly. Hydrodynamical aspects would 
be more important at the star-formation-scale, where e.g. 
fragmentation, winds and supernovae explosions play a cru- 
cial role. 

The galaxies were set on a prograde encounter with tidal 
tails developing on either side of the orbital plane, implying 
non-coplanar discfl The orbit of the associated point-mass 
problem is bound with eccentricity e ~ 0.96. The extended 
massive dark matter halos cause the galaxies to brake sig- 
nificantly already at the first pericentre passage (t ~ 11, 
snapshot used as our time reference). Tidal tails form soon 
after the first passage and continue to expand until the end 
of the simulat i on at t ~ 30. Comparing with HI maps of 
iHibbard et al.l (|200ll ) and optical images of the Antennae, 
the merger simulation gives a good match to the morphology 
and kinematics (peak systemic velocity of 130 km/s) at time 
t « 23, when a second passage through pericentre is in full 
swing (see Renaud et al. 2009, in preparation, for details). 



3 STATISTICAL APPROACH AND RESULTS 

Because the gas from which clusters form is largely confined 
to the disc of spirals, we focus on the fate of disc particles. 
The mass resolution of ~ 3 x 10^ Mq allows to treat indi- 
vidual disc particles as cluster-size element^. 

This resolution allows to monitor precisely the tidal field 
around a cluster. To study the effect of tides on a hypothet- 
ical star inside a cluster, let the cluster be at the origin 
of an inertial barycentric reference frame. The net force of 
the star is the sum of internal cluster and external galac- 
tic gravity. The tidal tensor computed from this external 
acceleration reads T'-' = dgl^t/dx-' , where x-' is the j-th 
component of the position vector and (;oxt is the acceler- 
ation. We retrieved the six-component tidal tensor of the 

4 X 10^ discs particles at each snapshot output time with 
a second-order finite difference scheme applied to a cube of 
side 6 = 40 pc. The symmetric tensor T admits a diagonal 
form with eigenvalues Ai ^ A2 ^ A3 computed using the 
method of lKopd (|2006l ). The signs of the eigenvalues deter- 
mine if the differential force leads to a stretch (A > 0) or 
a compression (A < 0) of the model cluster. Most examples 
of analytic potentials with a cored density show a mixture 
of bo th cases (Valluri 1993; Dckol et al. 2003; Gie les et al.l 
|2003). From here onwards, a fully compressive tidal mode 
is understood to be one where all three eigenvalues are neg- 
ative. In Fig. [1] we plot the percentage of disc particles 
in compressive mode as a function of time. The main evolu- 
tionary stages of the merger can easily be identified from the 

As the mass of the Antennae is measured with less accuracy 
than the velocities, we derived a mass-scale from = GM/r 
after adjusting the kinematics. 

^ See movies at http://astro.u-strasbg.fr/~renaud/simu.php 
^ Resolution requirements do not permit to simulate an individ- 
ual star clusters on the scale of a galaxy merger. 
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Figure 2. Left: map of gravitatio nal tides in the nuclear region of the simulated merger (~ 12 kpc wide). The black dots are candidate 
clusters identified bv iMengel et al.i (2005) from VLT images. The blue to red colours label extensive to progressively more compressive 
tidal modes. The two nuclei are marked by green crosses. Right: large-scale column-density map of the simulation at t ~ 23 time units 
(~ 300 Myr, peak 'D' in Fig.[TJ. 
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Figure 1. Percentage of the disc mass in compressive mode as a 
function of time. The peaks labelled 'A' through 'E' correspond 
to pericentre passages (A,D) and rapid morphological evolution 
(B,C,E), respectively. See text for details. 

rapid increase in particles undergoing compression. Peak 'A' 
occurs at the first pericentre passage, when the progenitors 
are largely whole, and the tidal field at its strongest. Around 
that time, bridges of matter form between the two progeni- 
tors and cause a second peak 'B' of compressive modes, of a 
somewhat longer duration. The number of compressive re- 
gions decreases steadily as the galaxies move apart, in the 
interval 50 < t < 250 Myr, prior to the second passage to 
pericentre. Once again a strong wave of compressive modes 



brings ~ 15% of the bodies to experience compressive tides 
(peak 'D' on Fig. [l|. Numerous short-duration spikes are 
visible around this second encounter (denoted 'C and 'E' 
on the figure) which are linked to streams of tidal clumps 
falling onto the newly formed core. 

Roughly speaking, the mass fraction in tidal compres- 
sion rises five fold when compared to the progenitor galax- 
ies in isolation. This relative increase compares well with 
hydrodynamical simulations of a major merger, where the 
number of star forming re gions also i ncreases on average five 
times at pericentre. iDi Matteo et al. I liooal report that 80% 
of a sample of ~ 900 GalMer merger simulations showed 
a relative star formation rate increase of ^ 5. Larger in- 
creases would require a specific configuration, or more ef- 
fective gas di ssipation (see e.g. iMihos fc Hernguistl Il996l : 
ISpringelll2000l ). 

Fig. [2] maps the spatial distribution of gravitational 
tidal modes of the nuclear region using a colour-coded 
scheme with blue meaning extensive tides, and red com- 
pressive ones. A projection of the particles on the sky is 
shown to the right for reference. The map was constructed 
with an 26.4 kpc x 26.4 kpc grid of 256^ mesh points. Large 
spiral patterns of compressive modes are clearly visible and 
reminiscent of the distribution of massive young clusters dis- 
patched throughout the tidal bridge linking th e two nuclei of 
the Antennae (the overlap region; see Fig. 1 of lMengel et al.l 
2005). An overlay of the tidal map to a grey-scale rendition 
of Mengel et al.'s image shows a clear correlation between 
the distribution of clusters (shown as dots) and the morphol- 
ogy and intensity of the compressive tides. The agreement 
with observations is especially good for the northern arc- 
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Figure 3. Distributions of the time-intervals At spent in com- 
pressive modes by disc particles, (a) top panel: differential frac- 
tion M{l\t + dt) — A1(At); (b) bottom panel: forward-integrated 
cumulative fraction A^(> At). The solid curves show the distri- 
butions of the longest uninterrupted sequences (lus), while the 
dash is for the total time (tt) in compressive mode. 



like feature and throughout the overlap region. To determine 
whether the coincidence of young clusters in regions of com- 
pressive tides is a chance event, we investigated the distribu- 
tion of their duration along individual orbits. Two statistical 
distributions were compiled, one for the Total Time {tt for 
short) spent in compressive tides during the whole simula- 
tion; the other covering the Longest Uninterrupted Sequence 
{lus for short) of snapshots with a particle immersed in a 
compressive tide. 

The results are graphed on Fig. [3] in log-normal axes. 
The top panel on the figure shows the differential mass frac- 
tion for both distributions, A4{At + dt) —A4{At), versus the 
time interval At; here dt = 2.5 Myr is the snapshot sam- 
pling rate. The bottom panel shows the forward-integrated 
distributions A4{> At), the mass fraction of particles ex- 
periencing a compressive tide of duration exceeding At, for 
each distribution. The four curves on Fig. |3] show a knee at 
At « 50 Myr {lus distributions) and 80 Myr {tt distribu- 
tions). Up to these times, the distributions are well-fitted 
with exponential decays A^(At) oc exp( — At/r) of charac- 
teristic times T ~ 8 Myr {lus distribution) and r ~ 20 Myr 
(tt distributions; cf. straight lines on Fig. [SJa). The interval 
r 10 Myr far exceeds the formation time of 0/B stars; 
what is more, about 40% of the lus sample experiences a 
compressive mode of duration longer than r = 10 Myr. The 
characteristic decay time r ~ 20 Myr of the tt sample com- 
fortably covers the formation phase of rich star clusters, in- 
cluding the residual-gas evacuation process. 

Because the model simulation discussed here involves 
gravitational dynamics only, a re-scaling leaving the dimen- 
sionless quantity GM/rv^ unchanged will give exactly the 
same results. Therefore, any major merger of mass unit AI' 
and length unit r' , will yield the same characteristic time- 
scale T ~ 10 Myr for the statistics of compressive tides pro- 
vided the ratio {M' /M)/{r' /rf = constant (cf. §2 for the 
units of Af and r). The same applies to the eigenvalues A of 
the tidal tensor. 



4 DISCUSSION 

The role of tidal fields in galaxy mergers is traditionally 
taken to be disruptive and to inhibit star formation. The 
analysis of a major merger simulation suggests instead that 
transient morphological features give rise to compressive 
modes of gravitational tides on scales ~ 40 pc, comparable 
e.g. to the cut-off radius of Milky Way clusters. The charac- 
teristic times r ~ 10 to 20 Myr derived for the log-normal 
distributions of duration interval At, Ai{At), serves to esti- 
mate the fraction of young Antennae clusters still embedded 
in compressive tides at the current time. When normalised 
so M{0) = 1, the quantity 1 - M{At = r) ~ 36% gives the 
fraction of clusters of ages ~ 10 Myr or less still confined to 
a region of compressive tide. These clusters would appear 
super- virialised which may explain the scatter seen in mass- 
to-light ratio among the younger population of Antennae 
clust ers (|Bastian fc GoodwinI l2006l : Ide Griis fc Parmentierl 
20011). Note that these clusters would all have formed during 
the second passage through pericentre (Fig. [1} , and the vast 
majority of them will soon emerge from their pressurised 
"cocoon" and experience the usual disruptive effect of the 
tidal field, causing a rapid increase in the dissolution rate of 
clusters (see the sharp drop around peak D, Fig. [3]). The 
statistics of compressive tides also imply that some two- 
thirds of clusters with ages > 30 Myr are now in extensive 
tides and may have suffered significant mass loss or total 
disruption: as a result the age distribution of young clus- 
ters should peak sharply ar ound ages ^10 M yr, causing a 
break in the power- law fit of lFall et all l|2005h . This conclu- 
sion is not inconsistent with their data for clusters in the 
mass range < 10^ M0. 

To contrast quantitatively the effect of tides on a star- 
forming region with the case of an isolated molecular cloud, 
we define the dimensionless ratio A = aXRy^ Rv/{GM tot) of 
tidal to gravitational binding energy of the cluster-forming 
cloud; A is once more the eigenvalue of a tidal tensor which 
is taken to be isotropic for simplicity, the virial radius 
of the cloud of total mass Mtot, and a ~ 1 a dimension- 
less quantity which depends on the mass distribution of the 
model cloud (e.g. a = 0.86 for a homogeneous sphere). Be- 
cause a fraction e < 1 of the total cloud mass is turned into 
stars, a young cluster forms stars while still embedded in 
residual gas. The ejection of this gas throug h feedback from 
stellar radiation may cause it to disrupt (iGoodwinI Il997l : 
iGe ver fc Burkertll200ll : iKroupa et al. ||2P01'). A classic result 
is tha t this will occur wheneve r e < 1/2 (Hills 1980; G oodwinI 
1997; Boi ly fc Kroupal l2003aV Repeating the derivation of 
iHillSi (jl980l l but with a background tidal field, we find a new 
equilibrium radius R = rR^ after gas expulsion, where r is 
obtained from 

2Ar3+(i-e-2A)r+|=0. (1) 

A solution r > still exists when e < 1/2 for all A < 0, 
corresponding to compressive tides. By contrast, a solution 
with A > requires e > 1/2 or more efficient star formation, 
as anticipated for a disruptive tidal mode. (The case A — 
gives the critical value e = 1/2 for isolated clusters.) Cru- 
cial to cluster survival, a compressive tide will hold together 
a star-forming region that would rapidly dissolve in isola- 
tion. This remains true even for a very low star-formation 
efficiency e. 
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iBoilv fc Kroupal l|2003bl ) performed N-body calcula- 
tions to address this question without any tidal field. Us- 
ing the same numerical set-up and an additi onal compres- 
sive tide with A = (a/5)A ~ -0.06, iFleckl IQ007) found 
a new equilibrium after expansion by a factor of ~ 2, for 
an SFE e « 0.10. Furthermore, the cluster would dissolve 
rapidly once the tide switches from being compressive, to the 
more likely extensive mode. This process therefore provides 
a unique clue to explain the discrepancy which exists be- 
tween the expected Gaussian and the near-power-law mass 
function in the Antennae. Th e weak dependance o f clusters 
half-light radius with mass (jvan den Berghl l2008l ) implies 
that compressive tides should be operative on all mass scales 
and severely imprint the mass function. As an aside, we note 
that the general character of compressive tides means that 
they are also at work at high-redshift, when merger rates 
were higher. 



5 CONCLUSION 

Using gravitational N-body simulation of the Antennae 
galaxy merger, we show that a significant fraction of the 
discs' mass undergoes totally compressive tides. This mode 
plays an important role by preventing a cluster from dis- 
solving. The compressive regions detected in the simulations 
show a very good match with observational data, both in 
term of their spatial distribution, and their characteristic 
duration which compare to observed clusters ages. 

Compressive modes of the tidal field help keep together 
clusters that would otherwise dissolve at formation, through 
e.g. rapid gas expulsion. By introducing a tidal term in the 
virial Eq. ([T}, we found that clusters should remain bound 
when they sit inside a compressive tide, even when the SFE 
is as low as 10%. The precise characteristics of the evolu- 
tion of clusters is clearly dependent on the relative strength 
of the tidal field. This is a hint that to understand the sur- 
vival of clusters, and especially the issues associated with the 
concept of "infant mortality", it is essential to include the 
full evolution of the environment in which clusters evolve. 
These super-virial young clusters sheltered by compressive 
tides will dissolve rapidly whenever the tides are no longer 
compressive. Hydrodynamical simulations will be required 
to link formally this effect with cluster formation histories. 
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